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Abstract 

I introduce a family of closeness functions between causal Lorentzian geometries 
of finite volume and arbitrary underlying topology. When points are randomly 
scattered in a Lorentzian manifold, with uniform density according to the volume 
element, some information on the topology and metric is encoded in the partial 
order that the causal structure induces among those points; one can then define 
closeness between Lorentzian geometries by comparing the sets of probabilities 
they give for obtaining the same posets. If the density of points is finite, one gets a 
pseudo-distance, which only compares the manifolds down to a finite volume scale, 
as illustrated here by a fully worked out example of two 2-dimensional manifolds 
of different topology; if the density is allowed to become infinite, a true distance 
can be defined on the space of all Lorentzian geometries. The introductory and 
concluding sections include some remarks on the motivation for this definition 
and its applications to quantum gravity. 
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I. Introduction 

The purpose of this paper is to propose a definition of closeness between Lorentzian 
geometries, where by Lorentzian geometry I mean a diffeomorphism equivalence class 
G = {(M, g)} of manifolds with Lorentzian metrics. More specifically, I will first define a 
pseudo-distance function d^{G,G') of two geometries G = {(M, g)} and G' = {(M', g')} 
with finite volumes Vj^^ and V^,, depending on an integer n, such that whenever d^{G, G') 
is small, the two geometries are close at large volume scales compared to Vj^/n and V^,/n, 
up to a global scale transformation; most of the paper is devoted to this pseudo-distance 
and its properties, but I will also extend the definition to a distance function d^{G,G') 
depending on a length parameter £. Notice that the geometries in question can be based 
on two entirely different manifolds M and M'. 

There are various contexts in which such a definition is useful, but the ones that mo- 
tivated this work are mostly related to quantum gravity. There is a growing amount of 
evidence, initially suggested by analogies with other theories and simple consistency argu- 
ments but increasingly supported by more rigorous results, that the structure of spacetime 
at the smallest scales (of the order of or smaller than the Planck volume £p — (Gh/c^)'^ 
— ^just for this equation, G stands for Newton's gravitational constant) differs significantly 
from that of the four-dimensional, topologically fiat differentiable manifold we use as a 
model in ordinary physics [1, 2, 3]. Very many different proposals exist for what to replace 
this manifold with; I will mention only a few of them here, as examples of situations in 
which one needs to talk about the closeness of Lorentzian geometries. 

If one assumes that large quantum fluctuations of the metric on small scales will be 
associated with fluctuations in the topology itself, but differentiable manifolds are still valid 
models for the geometry, one is led to the notion of spacetime foam [4,5,6], a bubbling 
topological magma in which topological entities like geons and wormholes fluctuate into 
and out of existence. Spacetime is a quantum superposition of differentiable manifolds of 
different topology, and the ones that contribute most to the classical spacetime we see are 
such that each topological fluctuation occupies on the average one Planck volume; at larger 
scales they are all thought to be close to each other, and essentially indistinguishable from 
a topologically flat manifold. 

On the other hand, there are hints that the very notion of manifolds and continu- 
ity may have to be abandoned for models that describe spacetime at Planck scales. In 
the causal set proposal, spacetime is considered as a locally finite partially ordered set 
[7, 8, 9, 10, 11, 12]; if the elements are thought of as events, occupying on the average one 
Planck volume each, the partial order is interpreted as giving the causal relations between 
them. In spin foam type proposals, the basic structure is also that of a graph, but with 
extra variables attached to the edges and vertices [13, 14, 15, 16, 17, 18]. In either case, the 
continuum and the rest of the Lorentzian manifold structure we see at large scales emerge 
as a thermodynamic limit, much like the description of a gas by thermodynamic quantities 
such as pressure and temperature emerges at large scales. Part of the reason why this 
limit exists is that, even though there are infinitely many Lorentzian manifolds which can 
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smoothly interpolate between the elements of a given discrete set, they are all supposedly 
indistinguishable at scales larger than £p. 

A definition of closeness between Lorentzian metrics on the same manifold M, in the 
form of a scale dependent function d;^(g,g'), has already been given in Ref 19 (the main 
idea can also be found in Ref 20). However, that definition is not diffeomorphism invari- 
ant, in the sense that, if (/> is a diffeomorphism of M, in general d^{g,g') ^ d-^{g, (j)*g'). 
In principle, given such a one can construct an invariant one [19], but in this case d-^ 
is difficult to work with; and, more importantly, it is not defined for metrics on different 
manifolds. My goal here is to set up a definition that is applicable to any two Lorentzian 
manifolds, analogously to the one given for Riemannian geometries by Gromov using geo- 
metrical concepts [21] or by Serin using spectral techniques [22]. Unfortunately, the ideas 
behind those distances rely heavily on the positive-definite nature of the metrics; the one 
I use here comes instead from causal set theory: G and G' are close if, when we distribute 
the same number of points at random with uniform density in (one representative of) each 
of them, the probability of obtaining any given induced partial order among those points 
is about the same in the two cases. A few of the ideas that led to this work appeared 
earlier in a different form in Ref 23. 

The use of uniform distributions of points is what makes the definition diffeomorphism- 
invariant, by not requiring us to identify points in the two manifolds; we are comparing 
instead the two geometries by independently sampling them, Montecarlo style, which brings 
a probabilistic aspect into the definition. Therefore, I begin in section II by briefly review- 
ing the definition and some properties of a uniform random distribution of points in a 
manifold, with respect to a given volume element. For simplicity, I will assume that all 
manifolds (M, g) have a finite total volume Vj^ = d^x \/—g-, where D is the dimension 
of M. If the manifolds have no closed timelike curves, each n-point sprinkling is endowed 
with a partial order by the causal structure on the manifold, and defines an element of 
the set of all partially ordered sets (posets) on n points. The idea then is to define 
d^{G, G') by comparing the two probabilities on corresponding to G and G' . The rest 
of the section is devoted to constructing a procedure for calculating those probabilities. 
Section III contains the definition of the family of pseudo-distances and a derivation of 
some of its properties, and section IV an example in which calculations are carried out in 
detail. In section V, I introduce a family of distances, which uses sprinklings of arbitrarily 
high numbers of points. The discussion is kept at a general level, independently of any 
applications, but the concluding section VI contains additional remarks on applications of 
this work as well as open issues. 

Finally, a few words concerning notation and terminology. Poset elements are denoted 
by P, q-i manifold points by x, y, ...; by the past or future of a point x in a Lorentzian 
manifold, I will mean its chronological past or future -f^(x) (this convention is adopted 
mainly for the sake of definiteness, since most of our considerations will depend just on the 
volume of those sets or of their intersections and unions, which for well-behaved geometries 
would be the same if I had used instead causal pasts/futures or their closures); and the 
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relationship y e /+(x), or x e I~{'y), will be indicated by x < y. Finally, P's will stand 
for probabilities, P's for probability densities, C's for posets, C's for sets of posets, Vj^ or 
V{R) for the volume of the region R C M, and R \ R' = {x. \ x E R, x ^ R'} for the 
difference between sets. 



II. Random Point Distributions and Partially Ordered Sets 

This section contains the elements that will go into the definition of the closeness measures. 
I begin with a summary of the few notions we will need regarding uniform distributions 
of points in a manifold, and then discuss how to obtain probabilities for different resulting 
partial orders. 

Given any manifold M with a volume element, in particular one with a metric (which 
at this point could be Riemannian or Lorentzian, possibly even degenerate — but not ev- 
erywhere, lest we get Vj^ = 0!), such that the total volume is finite, we can define a 
random process of sprinkling points uniformly by stating that, each time a point is chosen 
in M, the probability density that a particular x be picked is 

1 



PMi^\V^) = ^V-9i^), (1) 

in any coordinate system; equivalently, the probability that x fall in any given measurable 
region R C M (such as any interval or any finite union or intersection of such sets [24]) is 

Pm(x e P) = / Pm(x|x/=^) d"x = . (2) 

If the process is repeated n times, we get a imiform, random sprinkling of points with den- 
sity p := n/Vj^, or, if we forget the order they came in, an unlabelled n-point distribution; 
these are the events we are interested in, and for which we will calculate probabilities. 

One of the probabilities one uses most often in such cases is the one for exactly k 
points out of n to fall inside R (without specifying which ones) . This probability follows a 
binomial distribution, 

P{k, R I n, M) = (^^ ^^"^ ^ ^) ^ ^ \ 



which, as and n become very large, with p = const, approaches a Poisson distribution, 

P{k,R\n,M)^- ;f . (4) 

kl 

This last equation justifies the name Poisson distribution that is often used for the sets of 
points used in this paper, and corresponds to the infinite volume situation. The fact that 
in that case P{k, R \ n, M) can be written in the (exact) form (4), where only p appears 
and not n or V^, indicates that it may be possible to generalize the definitions and results 
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of this paper to infinite volumes, although in that case we do not have the probability 
density (1) available, which is what we would use to carry out an actual sprinkling, e.g., 
in a computer simulation. 

When we randomly sprinkle n points in M, the volume element ^/—g determines 
statistically where they will fall; given their positions, the causal structure g determines 
then the causal relations between them. Prom now on, all metrics will have Lorentzian 
signature and satisfy the past and future distinguishing condition (see, e.g., Refs 25 and 
26). In particular, this implies the causality condition (no closed causal curves), which 
guarantees that a partial order is induced on each sprinkling, defining an n-element poset 
C e C^; the slightly stronger distinguishing condition implies that "there are no almost 
closed causal curves," in a specific sense which gives some additional benefits, as I will 
discuss below. Diff'ercnt geometries G = {(M, g)} and G' = {(M', g')} will then in general 
give different probabilities P^{C\G) and P^{C\G') of obtaining each C E C^, which we may 
compare as a way to determine how close the geometries themselves are. It is therefore 
important to have a general procedure available for calculating, in principle at least, the 
probabilities P^{C\G). 

Let us start by fixing our notation. While is the set of unlabelled posets C on 
n elements, will denote the set of labelled n-element posets C, and E^(M) the set of 
n-point sprinklings a = (x^^, . . . , x^) in M. (One may argue that the labelling of the points 
should not be important; I am considering sprinklings to be ordered n-tuples of points here 
for convenience.) As already stated, our random events are n-point sprinklings a obtained 
as a result of a random process with uniform density. The volume element on M 

induces a probability density on E^(M); since the points are independently sprinkled, this 
can be obtained from products of single point probability densities (1) [27], 

P^(xi,...,x^|y=^) PM(Xilv^) = T7^rll-_i V-fK) ■ (5) 

If the spacetime (M, g) has no closed timelike curves, i.e., satisfies the chronology condition, 
the relation x-^ < X2 induced by the conformal structure g on M is a partial order, so the 
sprinkling a becomes a labelled poset C := | < pj iff x^ < x^- in cr}, i.e., we get a 
map $g : E^(M) given by cr C. This map is many-to-one, and the inverse image 

of any C is the set S — $r"^(C) C S^(M) of all sprinklings with the same induced labelled 
partial order. This set has non-zero measure in S^(M); in fact, its probability is 

^'e(^Iv^)= / PM(xi,...,xJv^)d-Xi...dX, (6) 

Js 

where S is specified by conditions on the relations between the sprinkled points giving, 
for each x^, a region C M it can fall into in order to have the right relations with the 
previously sprinkled x^ with j < i, according to C. Thus, the probability (6) is of the form 



This expression gives the probabihty that the sprinkhng give rise to a labelled poset C; we 
will see below how to specify the explictly. 

What we really want to find is the probability that the sprinkling give rise to an 
unlahelled poset C. Each C can be labelled in n\ ways, but in general some of these 
labellings are indistinguishable in terms of the order relation; more specifically, the number 
of permutations of elements of C that give the same C is the number of automorphisms of 
C, |Aut(C)| (this number is a property of C, independent of the specific C chosen), and 
we get that each C E can be obtained from n!/|Aut(C)| different labelled C's, so the 
probability we are looking for is 



|Aut(Oj| Vm -^-^^-l 7Mi(xi,...,Xi_i;C) 



(8) 



where C is an arbitrary labelling of C. 

Suppose a given labelling C = {p^} of C has been chosen to carry out the sprinkling. 
This means that, in order for the {x^} to be a realization of C, each must be in the 
future of the x^'s such that > pj, among the previously sprinkled ones, in the past of 
the ones such that < Pj, and spacelike related to the remaining ones. In other words, 
while X]^ can be anywhere, Mj^(C) = M, points x^ with i > 1 must fall in the regions 

_ _ _ f if P^ > Pj 

M,(xi,...,x,_i;C) = fl M,/x^.,C) , M,/x^.,C) = <^ /" if P, < P, (9) 

j<i M \Ij otherwise , 

where for futures and pasts I use the abbreviations if := /^(xj and I- := I~ U if . The 
most convenient labelling C to use in each case may vary. It is often a good choice to pick 
one compatible with the partial order on C, in the sense that if p- < p- then i < j, which 
can always be done (in fact, it just means "start labelling from the bottom and work your 
way up," and the choice is almost never unique); this has the advantage that, to reproduce 
the partial order on C, no x^ needs to be in the past of any of the previously sprinkled 
Xj 's with j < i, which eliminates the second case in (9). 

This completes the prescription for calculating the probabilities to be used in the 
closeness function. In practice, the dependence of each M- on the points x^, x^_j^ 
makes the probability very difficult to calculate analytically, and one would normally use 
other means such as computer methods, except for very simple situations like the one in 
section IV. 
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III. The Pseudo-Distance 

In this section, I will take the point of view that geometries G can be identified with the 
sets of probabilities {P^(C|G) | C e with a degree of approximation that improves as 
n increases. The task of defining a pseudo-distance between G and G' is then reduced to 
that of defining a distance between their respective sets of probabilities. I will do so, and 
then consider some properties of the resulting pseudo-distance. 

Various functions can be used as distances between sets of numbers; some simple ones 



to handle would be the ^^-type distance Sn\G,G') = \ Y.c^c^ l^n(C|<^) " ^n(C|<^')l> 



the Euclidean distance, or simply the "sup" distance, but in view of the interpretation 
of the numbers as probabilities, I will use instead the statistical distance introduced by 
Wootters [28] in the context of rays in Hilbert space, which is proportional to the number of 
statistically distinguishable, in an appropriate sense, intermediate probability sets between 
the two sets being compared. Let us then define, for any two geometries. 



as the coordinates of a point on the unit sphere, identifying a direction in probability space 
M'''"', and d^{G, G') is then proportional to the angle defined by the two corresponding 
directions; notice that, because all coordinates are non-negative, that angle is at most 7r/2, 
so with this definition d^{G, G') is at most equal to 1. 

Clearly, d^{G, G') is not positive-definite. For each n, the number \CJ\ of posets that 
can be made out of the n sprinkled points, although very large, is finite; thus, the value 
of d^{G,G') depends on a finite number of parameters, and cannot capture all of the 
information contained in the geometries. This means that d^ cannot be an actual distance 
function in the infinite-dimensional space of Lorentzian geometries. One possibility would 
be to take the limit n ^ oo; this may indeed give a distance, but it may be a trivial one, 
as I discuss below, and we shall consider a better alternative in section V. However, even 
for finite n, two geometries for which d^{G, G') — are close when probed at scales larger 
than the mean point spacing, and this is what we really need in some applications. 

Let us consider the other extreme situation, d^{G,G') = 1. For finite n, this can 
happen only for highly degenerate geometries, since it requires that the argument of the 
arccos function in (10) be zero, in other words that there be no C G for which both 
P^{C\G) and P^{C\G') are non-vanishing, i.e., which can be embedded in both geometries. 
One of the possible C's is always the totally ordered n-element poset (a chain), so one of 
the geometries (say, G) must assign zero probability to pairs of timelike related points; in 
G, the light cones of all points must have degenerated away to lines. Another possible poset 
is the totally disconnected one (an antichain), so one of the geometries (necessarily the 
other one, G') must assign zero probability to pairs of spacelike related points; G' has the 
wide open light cones of the infinite speed of light limit, or is a one- dimensional timelike 
line. No poset C e C„ other than those two can be embedded in either geometry. We 




(10) 



Geometrically, the fact that Y^^^q^ P^{C\G) = 1 means that y/P^(G\G) can be interpreted 
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conclude that the inequahty d^{G, G') < 1 cannot be saturated other than as a degenerate 
hmit of sequences of geometries of the type just described. 

In the hmit n oo, however, the situation may change. We know from contimmm 
results that the topology, differentiable, and conformal structures of a past and future 
distinguishing Lorentzian geometry can be recovered just from the knowledge of the causal 
relations between all pairs of points [29, 30] and that, if one considers instead pairs of 
points in a sequence of uniform sprinklings of increasing density, the same is true in the 
limit n — > oo, with the added bonus that the volume element can be recovered as well, up 
to a global factor [31, 20]; thus, in that limit, sequences of posets {C^}, where each C„ has 
n elements and is a subposet of the next one, C C'^+i, can be embedded at most in a 
single geometry G. This means that 

VGt^G", hm P„(CJG')PJCJG') = V{CJ. (11) 

n— »co 

In fact, it is also true that each individual probability P^{C^\G) or P^{C^\G') tends to 
zero as n ^ oo. But the number of terms in the summation in (10) grows very fast with 
n (faster than exponentially [32]), and the limiting value of d^{G, G') for G ^ G' depends 
on the rate of approach to zero of these probabilities. It is possible that d^{G^G') = 1 
for all G G' (in terms of the discussion above, many posets may be embeddable in both 
G and G' , but the limit is 1 because all products of probabilities in (11) go to zero fast 
enough) ; in this case would be a distance, but a trivial, not very useful one. 

We will see what d^{G, G') can be replaced by later in the paper, and return now to 
examing properties of d^ with finite n. In addition to its much greater ease of computation, 
the function d^{G,G') is also made interesting by the following reasonable conjectures: 
(i) In a sense, for large enough n, it is "almost a distance," or "positive-definite up to 
differences on small scales;" (ii) For any subset of geometries labelled by a finite number N 
of parameters (analogous to the "minisuperspaces" used for spatial geometries), there is a 
finite n such that d^ is a true distance function on this set, and (iii) For any two arbitrary 
(distinguishing, finite-volume) different geometries G and G' there is a finite n such that 
d^{G, G') > 0, with d^{G, G')^lasn^ oo. 

To start with, I will prove the intuitively obvious, and nice property of the closeness 
measure that it is a monotonically increasing function of n: 

\/G,G' dJG,G')<d^^,{G,G'). (12) 

To prove this inequality, consider the process of sprinkling n + 1 points in a geometry G 
as an n-point sprinkling, followed by the choice of one more point. Then, the probability 
of the first n points yielding any given C e C„ is a sum over probabilities for different 
C e obtained when the extra point is added, 

PJC\G) = P„,,(C'\G) , := ^ Q , (13) 
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where q, is the fraction of n-element subsets of C that are isomorphic to C, which 

can be expressed in terms of the number (^) of ways of picking an n-element subset of 
C that is isomorphic to C (this number may be called "C" choose C," and I will use the 
convention that it vanishes if C is not a subposet of C); notice that it is clear from the 
definition that YliC^Cr,, fc,C' = 1; for any C . Then, we can write 



2 

d {G,G') — — arccos 

TT 



(14) 

where it is understood that C E and C',C" G For each C, the corresponding 



term in the summation in (14) is of the form y (X^i ^j) where all and bj are 

non-negative, for which the general inequality 



holds. To prove this inequality, we can square the two sides, which gives a-bj and 

Si \A^A\A^7^' ^respectively; the terms with i = j are equal; separate the other ones 
in pairs, a-bj +(ijb- and 2-^ a^^ ^a^b^, respectively, and square them; since we always have 
af6| + 2a^ajb^bj + a^b^ > Aa^a^b^bj, (15) follows. Applying this to (14) gives 

dJG,GO < I arccos [E^ Ec^ ^{fc,C'?Pn+i{C'\G)P^^,{C'\G') 

= I [Ec (Ec J^c,c') ^JPn+i{C'\G) P„+i(C"|G') ■ 

= d^^,{G,G'). (16) 
As a consequence of the proof, we also see that 

d^{C,C') = d^+.iCC) iff d^+.iCC) = , (17) 

since the inequality in (16) can only be saturated if (15) is, and this will happen only if 
for all i and j, a^bj = a^b^, which in terms of our probabilities reduces to P^j^]^{C'\G) = 
P^_^i{C'\G'). As a byproduct, we also obtain the equality (13), which may be useful for 
calculating P^(C|G), or one of the P^_|_]^(C"|G)'s if the others are known. 
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IV. A Simple Example 

As an illustration of the definition of d^{G, G') and the procedure for calculating P^{C\G) 
introduced in section II, we consider a very simple example, which already involves two 
1-parameter families of geometries with different underlying manifolds: a finite-size rect- 
angular portion of 2-dimensional Minkowski space, with line element ds^ = —dt^ + dx^ 
and topology M ~ IR^, G^ = {(M, 77)}, and a similar one obtained after a spatial identifi- 
cation, with the same line element and topology M' ^ IR x S^, = {(M', ry)}. I will first 
introduce each geometry and calculate the simplest probabilites, ^2(^1^7) ^tnd ^2(^*1^^), 
then use these to find d2{G^,G'^); the results will give us an indication of features and 
limitations of ^2, and we will then see how to overcome these limitations by calculating 
the Pg's and using d^{G^, G'g). 

The geometry G^ is the rectangle M := {x | < x < a, < t < &} in two-dimensional 
Minkowski space. Since the probabilities we are looking for are invariant under a global 
rescaling, they cannot depend on the volume = a6, but only on the aspect ratio 
7 := b/a. For n = 2, C2 has two elements, the connected two-element poset [ and the 
disconnected one • •; we must calculate the integrals in (8) for these two posets. 

To get the connected poset J in a two-point sprinkling with the "bottom-up" labelling, 
we need X2 to fall in the future of Xj^, or M2(x;,^, C) = M2i(x;^, C) = and (8) becomes 

^2(1 \G^) = I d'x, / d^X2 = fdx fdt V{I+{x,t)) , (18) 

where x^^ = and the volume V{I~^{x,t)) is {b — t)^, with correction terms that are 

needed for some values of {x,t) (see Fig. 1), 

V{I+) = {b- tf - ~ ^ ~ e{x<b&xidt<b- x) 

_ {b + xf (19) 

(The step function d equals 1 if its argument is true, otherwise.) If we assume that 7 > 1, 
so that X <b and x > a — b are always satisfied, (18) gives 

1 -47 + 67^ 

To get the disconnected poset • • we need the points to be causally unrelated, M2(X]^, C) = 
M \ 7]^ , so (8) becomes 

P2{"\G,) = ^I d^x, / d% = l-P2(I|G) = ^^, (21) 

^ ■ J M Jm\Ii 07 

where instead of doing another integral I have used -P2(I 1^7) + -P2(* * [^7) — 1- 

For the case 6 < a, we can now either integrate (19) again, or use simple symmetry 
considerations. If we fiip the rectangle by exchanging a <-> 6, the manifold transforms 



^2(I|G,) = (20) 
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according to <-> Gj^^^; if we take the two sprinkled points -x.^ and X2 along, the posets 
are also turned into each other, Thus, 

2 (22) 
P,i' . \G^) = P,( I \G,/^) = ^ ^ . 

The geometry G'^ is the cylinder M' one obtains applying the spatial identification 
(t, 0) ~ {t, a) to the rectangle in G^, with the same line element; again, the probabilities 
only depend on the aspect ratio 7 := b/a. Similar calculations to the ones leading to (20) 
and (21), but now integrating 

V{I[+) ^ [h -tf - {h -t - \ af e{t <h - \ a) (23) 

(see Fig. 2) over M', give 



P^TI^M 1-67+127^ ( i^/N 67- 1 

^2(11^7) = ^2 ' P2{"\G^)^-^, (24) 



for 7 > |, and 

P2{\\G',) = il, P^{"\G'^) = l-li, (25) 

for 7 < |, when V{I['^) is just (6 — t)^; we cannot use a trick like the one in (22) here, but 
this probability is the easiest one to calculate anyway. 
If we now use the definition (10) to calculate 

d^{G^, G;) = ^ arccos P^iW ^) P^ii \G'^) + ^Jp^^^\gJpJ^^)] , (26) 

where the two geometries are characterized by the same parameter value 7, we get the 
function plotted in Fig. 3, which goes to zero as the aspect ratio 7 ^ or 00, and the 
difference between the two manifolds becomes immaterial because all pairs (X]^,X2) are 
spacelike or timelike related, respectively, in both geometries; but d2{G^,G!y) is not zero 
for any non-degenerate cases. However, if we use for geometries with different parameter 
values, we find that, for example, 

if7, <^<|, P,{l\G^)^P2{l\G's) when 5 = 7-172, (27) 

i.e., for each 7 there is a 5 such that d2{G^, G'§) vanishes; with two sprinkled points, the 
single available parameter P2(I|^) cannot distinguish all geometries in this example. We 
will now see that this can be done using three points. 

What we need to show is that, for all values of the parameters, among the elements 
ofCg = {•••,•!, V, A) {}) there is at least one which is embeddable in G and G' with 
different probabilities; because of property (12), we actually only need to do this for the 
parameter values for which d2{G^,G'g) = 0. Let us consider the three-element poset for 
which the probabilities are easiest to calculate, the linear order J (linear orders are always 
the easiest ones, because M^j{Xj,C) reduces just to and M^(xj^, ...,x^_^;C) to in 
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(9), if C is the "bottom-up" labelling). The calculations again proceed along the lines of 
those for (20) but are somewhat longer, since we now have to evaluate 

^3(! IG') - / d'x, / d^x^ / d2x3 . (28) 

^ - "^M Jm Ji+ Ji+ 

I will restrict myself to the case 7 < ^, where by explicitly integrating (28) one gets 

^'3(i|G,) = |7^-A7^ P3{l\G',) = h'. (29) 
When 5 = 7— 1 7^, a simple calculation gives 

P3il\G,)-Psil\G's) = jo^'i-3-h)^ (30) 
which does not vanish for 7 < |, so d^{G^,G'g) 7^ as expected. The two families of 
geometries are different enough that the induced order on a random three-element subset 
will pick out the difference. The other three-point probabilities P^{C\G) can be found 
without too much additional effort. I will now show how to do this, since the information 
those probabilities capture about the geometries is interesting in its own right, and because 
this will illustrate the use of some of the general relationships introduced above, as well as 
one new relationship and possible symmetries. 

Suppose we want to calculate the values of the probabilities P^{» • • \G), P^i'l \G), 
PsiY \G), PsiA \G), and P^il \G) for some geometry G, for which we already know the 
values oi P2{* *\G) and P2{1 \G); this includes having calculated V{I^) for every x^ e M, 
as in (19), or some other similar integral. It would be to our advantage to use as many 
relationships as possible among the P^{C\Gys. One is always given by the identity 



EceC3^3(C|G) = l, (31) 
and two more are always given by 

lPsi'l\G) + lPs{Y\G) + iP,{A \G) +Ps{l\G)=P2{l \G) , (32) 

P3(. . . |G) + I P3(. I |G) + i P3( V |G) + i P3( A \G) = • \G) , (33) 



arising from (13). Notice however that only two of the three relationships (31)-(33) are 
independent, since (13) already implies '^c'eCn 1 Pn+iiG'\G) = 1 if one uses valid P^'s and 
ffjQi^s^ satisfying X^cec Pn{G\G) = 1 and Xlcec fee ~ ^- additional relationship, 
not as simple as (32)-(33) but still useful, can be found by considering probabilities for 
subsets of C3 rather than just single P^{C\Gys. Since the outcomes V ^^d I are mutually 
exclusive, 

^'3(Vor 1|G') = P3(V|G') + P3(J|G), (34) 

which means that 

^^3(V|G)=P3(Vor \ \G)-P,{l\G) 

= 9^ I ^'^1 /.^'^^ /.^'^3 - ^^3(11^) 
^ ■ JM Jl+ Jl+ 

= 4- / {V{I+{x, t))f - P,{ [ \G) , (35) 
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where I have used the fact that, in a "bottom-up" labeUing, the condition for obtaining 
V or J is simply that both X2 and Xg be in the future of x^^. FinaUy, if G was time 
reversal invariant (the and G'^ in our example both are), we would get an additional 
relationship, 

P,{K\G)=P,{^\G). (36) 

We have found four relationships among the P3(C|G)'s; if they hold, only one probability 
needs to be calculated by direct application of (8). Also, in specific cases, it may be possible 
to use other symmetries of G to derive relationships of other types; for example, the one I 
used in the trick of (22) involving different parameter values. 

In our example, we start by calculating the integral in (35) using (19); if 7 < |, 

/ d^:>c{V{I+{x,t))f = r dx t dt{h-tf 
Jm Jo Jo 



+ 



+ 



[ dx [ 
Jo Jo 

[ dx f 

Ja-b Jo 



dt 



{h-t-xY 



b—a+x 



dt 



{b-t 



-{h-tf{h-t- xf 
a + x)^ 



-{h-tY{h-t-a + x) 



5 180 " 



(37) 



where one of the cross terms in the square of (19) does not contribute because in this case 
the intervals x < b and x > a — b don't overlap; then, using P^{1 \G) from (29), we find 
P^{Y\G), and from (36) we find P3( A 1^); with these, (32) gives P^{*l\G), and (33) gives 
P^{** * \G). Analogous calculations can be done for G', starting with 



/ d^x, {Vil+)f = r dx f dt {b -tf = \ ab' . 

JM' Jo JO 

To conclude, the full sets of probabilities we have calculated for 7 < |, are 



P2il\G,) = h-lr 



P2{"\G,) = l-ll+il' 
for two-point sprinklings, and 



P3i\\G,) = h'-T2^' 

Ps{V\G,) = h' 

P,{h\G,) = h' 
P3(.I|G^) = 27- 

i^3(---|G,) = l-27+S7^-M7 



^2(1 1^;) = 

P2{"\G'^) 
P,{\\G')- 



37 
1 - 



7 



- 1 ^3 
5 ' 

- 1^3 



10 



5 ' 



(38) 



17 ^,3 



1^2 

^^3(71^;) = §7^ 

P,{h\G',) = li^ 
P3(.I|G;)=27- 
P3(...|G') = l-27+f72 



for three-point sprinklings. The corresponding ones for 7 > | can be similarly calculated 
with the formalism described above. 
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V. The Distance Function 

We can now extend the definition of the closeness function to a distance. Only the definition 
and a few comments will be given here; a more extensive study of its properties is left for 
future work. From the previous discussion, it should be clear that in this case we need 
to let the number of points sprinkled in each manifold go to infinity, so that we probe its 
structure at arbitrarily small scales. Also, it is not sufficient to let n — oo in d^{G,G'), 
both because the resulting distance may be trivial, and because, like all (i^'s, it would not 
distinguish between different values of the total volumes Vj^ and V^, . Each random event 
considered in previous sections was the choice of n points in each manifold; since n was 
the same for both manifolds, the outcomes gave us no information on their total volumes. 
We can overcome this limitation in the distance by letting each random event consist in 
the choice of both n and the location of the points, and this will give me an opportunity to 
mention one feature of the closeness functions that had not explicitly come up until now. 

Consider two geometries G = {{M,g)} and G' = {{M' , g')}, as before. We will draw 
the number of points n to be sprinkled in each manifold from Poisson distributions, whose 
means will be proportional to the respective volumes, and will thus in general be different, 
the relationship between the two means being given by the point densities they correspond 
to. However, the two manifolds may be of different dimensions, D and D' respectively, 
in which case it would be meaningless to require the two volume densities of points to be 
equal; what we can require is equality of the "mean point spacings", i.e., that the volume 
densities p and p' satisfy (p)^^^ = {p')^^^ = (For example, in quantum gravity 

applications, wc can think of £ as being the Planck length £p, and the issue of different 
dimensionalities is relevant for higher-dimensional theories such as the Kaluza-Klein ones 
— for a recent review, see Ref 33 — where we might want to compare a macroscopic four- 
dimensional manifold to a D-dimensional fundamental one which is, at least locally, 
considered to be a product of the type W ~ %f x ^"^M, with ^"^M of volume £^~'^.) 

To define the distance function, choose a positive mean point spacing i around which 
most of the contribution to the distance will come from; sprinkle points in G and G' by 
first choosing, each time, the number of points according to Poisson distributions 

respectively, where := Vj^/£° and p' :— Vj^,/£°' , and distribute in {M,g) and {M',g') 
the chosen numbers of points uniformly at random; the probabilities of obtaining any given 
poset C e result are now, respectively, 

P,{n,C\G) = P^{n)P^{C \ G) , Pe{n,C\G') = P^,{n)P^{C \ G') ; (40) 
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finally, compare these probabilities by extending (10) to 

oo 

E \/Pein,C\G) ^P,{n,C\G') 

n=0 C€C„ 



2 

d^{G,G') : = — arccos 

TT 



2 

— arccos 

TT 



(41) 

Ln=0 ^ cec„ ' ^ . 

Here, I am assuming we have defined the probabilities P^{C\G) for the sets of one-element 
posets, Ci = {•}, and zero-element posets, Cq = {0}; if we set Pi{*\G) — 1, consistently 
with the general definition, and adopt the convention that Pq{0\G) = 1, the argument of 
the arccos function can be written as 

^e-Me-M' + y'/i/i'e-Me-^' + 5^(^P^(n)P^,(n) ^PJC|G) PJC|G') ) • (42) 

The expressions (41) and (42) are clearly well-defined; the rapid decrease of P^in) and 
P^i (n) for large n makes them finite, and the fact that they are probabilities implies that 
(41) actually gives a number between and 1, for the same reason as the pseudo-distance 
(10) did. To examine these two extreme situations, consider (42). This expression vanishes 
only in the large ^ or ^' limit (£ 0, or at least one of the volumes oo), so that the first 
two terms vanish, and if the contribution from all n > 2 vanishes; this may imply that the 
conformal structures have the degeneracies described earlier in section III for d^{G,G'), 
at least if the manifolds have equal volumes, but in any case we can already see that 
d(^{G, G) = 1 only in situations obtained as limits of ones of the type under consideration. 

The more interesting situation is when d^{G, G') = 0. We can see from (41) that this 
implies P^{C\G) = P^{C\G') for all n and C, since for all n the summation over C E 
must equal 1, and P^{n) = P^,{n) for all n, since the sum over n also must give 1. The 
latter conditions obviously mean that Vj^ — Vj^,, and the former set implies G — G'; the 
sketch of a proof goes as follows. When we sprinkle an increasing number of points in a 
manifold, we build a sequence of posets {C'^j^eN' with C C„_|_i for all n. In the limit 
n — > oo, we obtain the direct limit of this sequence, := {J^^^C^, and the equality 
of all P^{C\G) = P^{C\G') turns into the equality of appropriately defined probabilities 
P^{C\G) and P^{C\G') on C^. Now, if we apply a completion procedure to C^, analogous 
to the Dedekind cut construction of real numbers from rational ones, we get the points of 
the original manifold together with their causal relations and the conformal factor, i.e., we 
get back the geometry G [20, 31]. But if the two probabilities P^{C\G) and P^{C\G') are 
equal, and the completion of infinite posets drawn from them gives respectively G and G', 
with probability one, it must be that G = G', and thus di is positive-definite. 
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VI. Concluding Remarks 

To measure the closeness of Lorentzian geometries, I have introduced a family of pseudo- 
distances d^{G, G') and a family of distances G') on the space of all past and future 
distinguishing Lorentzian geometries of finite volume. The main idea was to sprinkle 
points uniformly at random in G and G', and use the resulting probabilities P^{C\G) as 
the basic ingredients for the functions; in this paper, those probabilities were combined 
using one specific distance between probability measures, but others are known and may 
be more suitable in some applications. The closeness functions presented here, together 
with other possible such functions based on the same probabilities, are the only non-trivial 
diffeomorphism-invariant ones on this space that I am aware of. 

A number of interesting questions arise about the statistical approach to Lorentzian 
geometry discussed here. Even before their use in the closeness functions, the probabilities 
P^{C\G) are interesting in themselves, as a complete set of invariants (together with the 
volume) of finite volume, distinguishing Lorentzian geometries. It would be worth while to 
study the type of information about the manifold that those different invariants contain; 
for example, how they encode dimensionality, how they are affected by conformal transfor- 
mations as opposed to changes in the conformal structure, or how one can tell "localized" 
changes from "global" changes in a manifold from their effect on the P^{C\G). Possible 
starting points in answering these questions may consist in examining examples along the 
lines of the one in section IV but in which different parameters are varied, e.g., comparing 
a two-dimensional and a three-dimensional manifold, or modifying one by a conformal 
transformation; and studying analytically the effect of small variations g i-^ g + 5g. 

The answer to questions of the above type may then allow us to word in a more 
precise way statements like "geometries for which d^{G,G') is small are close down to 
the scale Vj^/nf understand how the topology induced by the ci^'s and d^'s on the space 
of Lorentzian geometries relates to previously studied ones [34] ; and place bounds on the 
value of d^{G, G') when the actual value cannot be calculated, for example through bounds 
on the probabilities due to the non-embeddability of some C's in a geometry, such as 
C"s that require higher dimensions. The infinite density limit and the properties of d^^ need 
to be understood better than what is sketched in section V, and the limit in which the 
"regulator" (. is taken to zero is a potentially useful one. It would also be useful to extend 
the present work to a definition of closeness that applies to infinite volume manifolds, as 
mentioned in section II; in that case, one may need to introduce a quasi-local element in 
the definition, and use finite size subsets of sprinklings of density p. 

On the physical side, this work may be related to definitions of approximate solu- 
tions of Einstein's equation [35, 36, 37], and spacetimes with approximate symmetries [38], 
which have been considered for various reasons, including their relevance to the issue of 
gravitational entropy and the smoothing problem in cosmology [39]. These problems, in 
addition to the motivation coming from quantum gravity, make it an interesting issue to 
study properties of P^(C|G), by analytical methods or numerical simulations. 
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Figure Captions 



Figure 1: The geometry G^. The drawing shows the case 6 < a, or 7 < 1, with a 
sprinkled point x^^ and its future hght cone. For this particular point, y(/+(X]^)) = 
(6 — t)"^ — \ [{b — t) — [a + ic)]^, where the area of the small triangle with dashed 
edges in the upper right hand corner must be subtracted, since a — x < b — t. 

Figure 2: The geometry G'^. The drawing shows the case b > a/2, or ■y > ^, with a 
sprinkled point Xj^ and its future light cone. For this particular point, K(/"'"(xj)) = 
{b — t)'^ — {b — t — a/2)^, where the area of the two small triangles with dashed edges 
must be subtracted, since b — t > a/2. Any two points whose x coordinates differ by 
a are identified; in particular, the two outer vertical dashed lines are to be identified 
with each other. 

Figure 3: Plot of d2{G^,G'^) as a function of 7, for < 7 < 5; for 7 > 5, the 
function decreases monotonically, and approaches zero as 7 — > 00. 
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